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ABSTRACT 

Microwave Sounding Unit (MSU) Ch 2 data sets, collected from sequential, polar-orbiting. Sun- 
synchronous National Oceanic and Atmospheric Administration (NOAA) operational satellites, 
contain systematic calibration errors that are coupled to the diurnal temperature cycle over the 
globe. Since these coupled errors in MSU data differ between successive satellites, it is neces- 
sary to make compensatory adjustments to these multisatellite data sets in order to determine 
long-term global temperature change. With the aid of the observations during overlapping 
periods of successive satellites, we can determine such adjustments and use them to account for 
the coupled errors in the long-term time series of MSU Ch 2 global temperature. In turn, these 
adjusted MSU Ch 2 data sets can be used to yield global temperature trend. 


In a pioneering study, Spencer and Christy (SC) (1990) developed a procedure to derive the 
global temperature trend from MSU Ch 2 data. In their procedure, the magnitude of the coupled 
errors is not determined explicitly. Furthermore, based on some assumptions, these coupled 
errors are eliminated in three separate steps. Such a procedure can leave unaccounted residual 
errors in the time series of the temperature anomalies deduced by SC, which could lead to a 
spurious long-term temperature trend derived from their analysis. In the present study, we have 
developed a method that avoids the shortcomings of the SC procedure. Based on our analysis, we 
find there is a global warming of 0.23 ±0.12 K between 1980 and 1991. Also, in this study, the 
time series of global temperature anomalies constructed by removing the global mean annual 
temperature cycle compares favorably with a similar time series obtained from conventional 
observations of temperature. 
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1. INTRODUCTION 

Temperature observations made from satellites orbiting 
around the Earth can be valuable to studies of global 
temperature and its variation with time, because of the 
uniformity of global measurements made by them. How- 
ever, necessary accuracy in the data and in the method of its 
analysis should be maintained to get good results. With 
such measurements, global temperature could be monitored 
readily, which could help assess the impact of anthropo- 
genic effects. 

Passive infrared radiometer temperature observations made 
from satellites are contaminated by clouds in the tropo- 
sphere. Therefore, these observations are not helpful to 
monitor global temperature change, which can be on the 
order of 0.1 K in a decade. In a pioneering study, Spencer 
and Christy (SC) (1990) demonstrated in a preliminary way 
the potential of the passive microwave radiometers to 
monitor global temperature with that precision. In particu- 
lar, they emphasize that Microwave Sounding Unit 1 Ch 2 
near 53.74 GHz, because of its weak sensitivity to hydro- 
meteors in the atmosphere and very strong response to 
thermal state, provides measurements that reflect basically 
the temperature of the midtroposphere, which is closely 
related to the surface air temperature. 

The global temperature trend deduced from MSU data by 
SC disagrees with that obtained from conventional data 
analysis, which is mainly based on surface air temperature 
or radiosonde-observed temperature (Hansen and Lebedeff, 
1987; Jones et al., 1986; Angell, 1988; Oort and Liu, 1993). 
The global temperature trend derived by SC (1992) from 
MSU Ch 2 data is near zero (-0.03 K) during the time 
period 1979-1990. On the other hand, the conventional data 
analyses show a wanning of about 0.2 K during this period. 
This disagreement has led to a careful examination of the 
analysis of SC by Hansen and Wilson (1993), Trenberth and 
Hurrell (1996), Prabhakara et al. (1995), and Shah and Rind 
(1995). Utilizing the radiosonde data, Hansen and Wilson 
simulated theoretically a time series of Ch 2 global tem- 
perature from 1979 to 1991 and calculated the resulting 
temperature trend. Based on their study, they indicate the 
analysis of SC underestimates the global temperature trend. 
Trenberth and Hurrell contend that the limited amount of 
overlap data used to calculate intersatellite calibration 
corrections can introduce errors in the analysis of SC. 
Prabhakara et al. suggest that there is nonnegligible residual 


hydrometeor contamination in the MSU Ch 2 data, which 
has not been screened by SC; and that this residual contami- 
nation could effect the resulting trend. Shah and Rind 
emphasize noise introduced by the surface conditions. The 
above critical examinations did not reveal the key weakness 
in SC's satellite data analysis. In this study, we find there is 
a weakness in the SC treatment of the diurnal effects and 
the calibration errors. When this weakness is removed, the 
satellite observations indicate a global warming in the 
period 1980-1991 of about 0.23 ±0.12 K. 

If the long-time series of successive NO A A operational 
satellites had the same Sun-synchronous orbital geometry 
over the globe and contained no errors, it would be a 
relatively simple matter to get the time series of the global 
mean temperature from MSU Ch 2 data. However, this is 
not so. There are two different orbital geometries for 
NO A A satellites. One of them has local equatorial crossing 
times set close to 7:30 a.m. and 7:30 p.m., while the other 
has these local times set close to 2:30 a.m. and 2:30 p.m. A 
satellite with the first orbital geometry is referred to as a 
“morning satellite,” while a satellite with the second 
geometry is called an “afternoon satellite.” 

Sometimes, a morning satellite is followed by an afternoon 
satellite. This introduces significant problems in developing 
a time series of MSU Ch 2 temperature data that can be 
analyzed to get a long-term trend. First, there can be 
differences in the Ch 2 temperature of the successive 
satellites because of the diurnal cycle. Secondly, the 
calibration of the MSU data between the two successive 
satellites could differ. These two errors exist in the MSU 
data as coupled errors. To eliminate these coupled errors , 
Spencer et al. (1991) developed a procedure which depends 
crucially on the temperature anomalies with respect to a 
suitable mean annual cycle in Ch 2 temperature data. The 
three steps involved in this procedure are described below. 

First, for each latitude belt the NOAA 6 1982 mean annual 
cycle is subtracted from the MSU Ch 2 data of all the 
morning satellites to get the zonal temperature anomalies of 
the morning satellites, while the NOAA 7 1982 mean 
annual cycle is subtracted from the MSU Ch 2 data of the 
afternoon satellites to get the temperature anomalies of the 


1 Microwave Sounding Unit flown on Sun-synchronous, polar- 
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afternoon satellites. These anomalies are properly weighted 
and integrated to get global anomalies. It is assumed that 
the diurnal cycle effect is removed in this step. In the 
second step, the temperature anomalies obtained from the 
overlapping data from two successive satellites are equal- 
ized with the help of an adjustment term. From this step, it 
is assumed that the intersatellite calibration errors are 
removed. In the final step, the temperature anomaly time 
series obtained as a result of the two previous steps is 
refined by subtracting an improved mean annual cycle. This 
procedure involves a chain of subtractions to get the final 
anomaly time series. We find that the long-term global 
temperature trend derived from this final time series does 
not account fully for the coupled errors in the MSU Ch 2 data 

The objective of the present study is to deduce global 
temperature from the MSU Ch 2 data, using an independent 
approach. For this purpose, we have to understand the 
nature of the coupled errors present in the MSU data, so 
that we can explicitly adjust the data for these errors. Then 
we will develop a method, based on the MSU Ch 2 global 
temperature data (not the temperature anomalies ) , to 
determine the trend. 

2. DATA ANALYSIS 

2.1 Selection of data 

The MSU is a cross-track scanning instrument that makes 
11 observations along each scan (NOAA, 1991). In the 
present study, we have selected only the measurements 
made by the MSU in the nadir direction. This choice 
reduces the number of observations, but it also avoids the 
need to empirically relate the cross-track scan data to the 
data along nadir direction, as done by SC (1992). 

Outliers in the MSU Ch 2 brightness temperature data are 
eliminated. Whenever this temperature exceeds 300 K, or is 
less than 200 K, we neglect such data. Such outliers 
constitute less than 1 % of the entire data set. 

For each month from 1980 to 1991, the MSU nadir data has 
been binned into 2° latitude x 3° longitude grid boxes over 
the globe. The data set from the morning satellites is 
separated into two global sets by the satellite equatorial 
crossing local time. Thus, one set corresponds closely to the 
7:30 a.m. equatorial crossing local time and the other set 
closely to the 7:30 p.m. equatorial crossing local time. 
Similarly, the afternoon satellite data set is partitioned into 


two global data sets by the 2:30 a.m. and the 2:30 p.m. 
equatorial crossing local times. This step of the analysis 
reveals that, because of the orbital configuration of the 
polar-orbiting NOAA operational satellites, for each 2° x 3° 
grid box near the equator there is a maximum possible 
number of 7 MSU Ch 2 nadir observations in a given month 
at each of the four local equatorial crossing times. This 
number increases away from the equator. Subsequently, a 
monthly average of these observations, for each month and 
in each grid box, is formed separately for the 7:30 a.m. and 
the 7:30 p.m. local times for the morning satellites. Simi- 
larly, for the afternoon satellites, grid box monthly averages 
are formed separately for the 2:30 a.m. and the 2:30 p.m. 
local times. 

For the purpose of our analysis, the MSU data in the 2° x 3° 
grid boxes are separated into global land and ocean areas. 
Data from the grid boxes that form the coastal regions are 
not considered in the analysis. The data loss due to coastal 
rejection does not exceed 20% in the 60N to 60S zone, as 
shown in figure 1 ; and only in limited parts of the high 
latitudes does the rejected coastal data exceed 20%. The 
total land plus ocean area considered in our analysis 
amounts to 84% of the area from 75N to 75S, of which the 
land accounts for 27% and the ocean accounts for 57%. 

This rejection of coastal data does not bias the global 
proportion of ocean and land areas. 

We have not eliminated any of the MSU Ch 2 data for 
hydrometeor contamination, because the interannual error 
in MSU Ch 2 data introduced by this contamination over 
the globe, -0.04 K (Prabhakara et al„ 1995) is much smaller 
than other errors in the MSU data. 

There are deficiencies in the Ch 2 global monthly mean 
temperatures of the following months: March 1980, June 
1982, Dec. 1982, Sept. 1984, Feb. 1988, Dec. 1988, and 
Feb. 1989. These deficiencies arise due to unusually long 
time gaps in the MSU data set, or when there are significant 
numbers of questionable data within a given month that are 
close to, but within, the limits of the outlier rejection. Based 
on first and second order derivatives, these deficiencies are 
remedied by interpolation of the time series of the Ch 2 
global monthly mean temperature. 

The monthly population of the Ch 2 nadir data in a grid box 
for a particular time of the satellite orbit may not appear 
large; but we find it is adequate. This is due to the fact that 
our analysis depends on multiyear average temperature, and 
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not the anomaly of the monthly average temperature, on a 
global scale. Temperature change over a long period, or its 
time trend, is deduced in our study from the multiyear 
averages over the global ocean and the global land. 


superimposed on the calibration errors, at present the 
magnitude of either the diurnal effect or the calibration 
errors cannot be assessed individually. For this reason, 
these coupled errors have to be treated as one entity. 



Figure 1 . Percentage of Land, Ocean, and Land-Ocean 
Sum in each latitudinal belt from 75N to 75S. The percent- 
age in each belt constituted by coastal grid boxes is given 
by (100 - Land-Ocean Sum). 

2.2 Nature of the coupled errors in MSU data 

A thorough discussion of the MSU calibration procedure, 
utilized to get NOAA operational data sets, is outside the 
scope of this study, but a brief description of this procedure 
is given here. In this calibration procedure, the temperatures 
of the onboard warm black body (-288 K) and free space 
(-2.7 K) are first taken as references. Then, a given radiom- 
eter reading of a target on the Earth is converted into a 
brightness temperature using a calibration procedure 
provided by NOAA (1997). 

From the earlier description of the satellite orbital geometry 
of the morning and afternoon satellites, we note that 
whenever both morning and afternoon satellites are observ- 
ing the Earth, such observations made by those satellites in 
a month can be partitioned into four sets based on each 
local equatorial crossing time of the orbits. We expect that 
these four sets of data will reflect the phase and amplitude 
of the global diurnal temperature cycle, if there are no 
calibration errors in the observations of the two satellites. 
However, since there are calibration errors, the diurnal 
temperature cycle in these data will be distorted. Because 
the variations of temperature in the diurnal cycle are 


In order to understand the nature of these coupled errors , 
we have analyzed the MSU Ch 1 and Ch 2 temperature 
data, obtained from overlapping NOAA 6 and NOAA 7 
satellites, separately over the global land and ocean areas 
for 1982. MSU Ch 1 observations (50.3 GHz) lie in the 
window region of the 60 GHz oxygen absorption band, so 
they can provide information that supports and supplements 
the information given by Ch 2 data. The MSU Ch 1 or Ch 2 
global monthly mean temperature T^'for a given period is 
calculated using the following numerical equation: 

, X i( T AM ) ij W iJ C ° S ( 2 l i l) 

Tm< jsl . (I) 

37 120 

X 5> SJ cos(2|j|) 

i=— 3 7 j=l 

In equation 1 , (T^). is the monthly mean Ch 1 or Ch 2 
temperature at the grid box (i, j) measured close to 7:30 
a.m. or 2:30 a.m. The index i, when multiplied by a factor of 
two, represents the central latitude of the 2° latitude x 3° 
longitude grid box. As i is incremented from -37 to 37, this 
latitude increases from -74 to 74 in steps of 2°. The grid box 
that has the equator as its central latitude is given by i = 0. 
The index j represents the central longitude of the grid box. 
As j increases from 1 to 120, the longitude increases in 
steps of 3° from 178.5° W to 178.5° E. Weighting by the 
cosine of latitude is done to get proper area average over the 
spherical Earth. The weight w assigned to each grid box is 
normally equal to one. However, when we are estimating 
the global Ch 1 or Ch 2 mean temperature over lands, w y is 
set equal to zero for the grid boxes representing oceans and 
coastlines. Similarly, when we are estimating Ch 1 or Ch 2 
mean temperature over oceans, w„ is set equal to zero over 
lands and coastlines. T m 'is calculated in a similar fashion 
from the data observed close to 2:30 p.m. and 7:30 p.m. 
local time. In this study, we are interested in the average of 
T^'or Trover global land or ocean for time periods that are 
much longer than 1 month. Such an average for K number 
of months is given by 


Tam 



( 2 ) 
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As stated previously, when we separate the 1982 overlap- 
ping MSU Ch 1 observations made by NOAA 6 and NOAA 
7 into four sets based on the local equatorial crossing time, 
we can display graphically the annual temperature cycle 
over global land or ocean at four local times of the day. In 
figure 2a, such annual temperature cycles over the land are 
shown. In Ch 1, there is strong information on the surface 
temperature and the surface emissivity. From the four 


curves shown in this figure, one can get a rough estimate of 
the amplitude and phase of the diurnal temperature cycle for 
each month of the year, mainly because the amplitude of the 
diurnal cycle on the land, -10 K, is large compared to the 
calibration errors. In figures 2b and 2c, we show the diurnal 
temperature cycle in January 1982 and July 1982 based on 
the temperature at the four local times of the day on the 
global land. 


MSU Ch 1 Annual Temperature (T b1 ) Cycle 

Over Global-Land for 1982 
Derived From NOAA6 and NOAA7 
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Figures 2a, 2b, and 2c, Annual temperature cycle over the global land for 1982 deduced from 
NOAA 6 in MSU Ch 1 at 7:30 a.m. and 7:30 p.m., and from NOAA 7 for the local times 2:30 a.m. 
and 2:30 p.m. b) Diurnal cycle over global land for January 1982, and c) for July 1982. 
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Similar graphical display of Ch 1 data over global ocean is 
presented in figures 3a, 3b, and 3c. Note that because the 
sea surface emissivity in the microwave region is ~0.5, the 
brightness temperature of Ch 1 over the ocean is much 
smaller than that over the land. From these figures, it can 
also be seen that the amplitude of the diurnal cycle over 
ocean is small. Furthermore, we see that the MSU Ch 1 data 
indicates a minimum near 2 p.m., which is contrary to the 
Tropical Ocean Global Atmosphere-Coupled Ocean 
Atmosphere Response Experiment (TOGA-COARE) 


observations of the diurnal cycle over the equatorial Pacific 
region. TOGA-COARE observations show a maximum in 
the surface temperature close to 2 p.m. and the diurnal 
amplitude is about 1 K (Sui et al., 1997). Thus, we infer the 
calibration errors in Ch 1, in particular those near 2:30 p.m., 
are large enough (-0.5 K) to badly distort the weak diurnal 
cycle present over the oceans. Much smaller systematic 
calibration errors that could be present in Ch 1 data at other 
local times of the day are not readily noticed in figures 3b 
and 3c. 


MSU Ch 1 Annual Temperature (T b1 ) Cycle 

Over Global-Ocean for 1982 
Derived From NOAA6 and NOAA7 



26.5 

January Diurnal Cycle 
Over Ocean 

7 

July Diurnal Cycle 
Over Ocean 

26 

1 1 1 » 

: (b) 

~ 26.5 

; — o 1 1 1 — ; 

: (O o 

g 

8 25 5 

• 0 o - 

r o - 

I 263 

o 

L 25 

T o' 

H S 26.1 

o 

H 

24.5 

L — i i i i : 

25.9 

1 1 1 1 


230am 730am 230pm 730pm 

230am 730am 230pm 730pm 


Time of Day Time of Day 


Figures 3a, 3b, and 3c. Same as figures 2a, 2b, and 2c, except for MSU Ch 1 data over global ocean. 
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In figures 4a and 5a, we present a graphical analysis of 
MSU Ch 2 data on global land and ocean similar to that 
presented in figures 2a and 3a. The diurnal temperature 
cycle in Ch 2 over land and ocean in January 1982 and July 
1982 are presented in figures 4b and 4c, and figures 5b and 
5c, respectively. As in the Ch 1 data, the characteristics of 
the calibration errors can be found in the Ch 2 data. We 
point out that the amplitude of the diurnal cycle on the lands 


in Ch 2, shown in figure 4a, is much weaker compared to 
that in Ch 1. The erroneous 2:30 p.m. minimum, similar to 
that in Ch 1, is also noticed in the diurnal cycle of Ch 2 over 
ocean (figures 5b and 5c). Explanation of all these system- 
atic errors is beyond the scope of this study. However, we 
believe these errors could be related to instrument exposure 
to sunlight. Since these errors are systematic, it is possible 
to make suitable compensatory adjustments for them. 


MSU Ch 2 Annual Temperature (T b2 ) Cycle 
Over Global-Land for 1982 
Derived From NOAA6 and NOAA7 
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Figures 4a, 4b, and 4c. Same as figures 2a, 2b, and 2c, except for MSU Ch 2 data over global land. 
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MSU Ch 2 Annual Temperature (T b2 ) Cycle 

Over Global-Ocean for 1982 
Derived From NOAA6 and NOAA7 
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Figures 5a, 5b, and 5c. Same as figures 2a, 2b, and 2c, except for MSU Ch 2 data over global ocean. 


Analysis of MSU Ch 1 and Ch 2 data, and the diurnal cycle 
presented above, could be useful also in future studies. 
Utilizing radiative transfer theory, we can decouple objec- 
tively the calibration errors from the diurnal cycle with the 
help of radiosonde observations. Such radiosonde data 
should correspond to the time of satellite observation and 
should be representative of the large footprint (-100 km) of 


the radiometer. Also one should know the surface tempera- 
ture and microwave surface emissivity. At this time, we do 
not have such measurements of the atmosphere and the 
surface. As an alternative, we seek a crude solution to this 
problem, which is elaborated in the next section. In this 
alternative solution, the coupled errors are eliminated 
jointly over an area of global scale. 
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2.3 Rectification of MSU Ch 2 data 

Basically, in the MSU Ch 2 data, the annual cycle and the 
diurnal cycle of temperature are intertwined. These cyclical 
effects need to be removed. In this study, elimination of the 
annual cycle is treated in a simple way. This is done by 
averaging the data over two or three complete annual 
cycles. The MSU data up to the end of 1991 , obtained from 
NOAA 6 to NOAA 11, allows us to define such multiyear 
averages for each satellite starting from 1 January 1980 and 
ending with 31 December 1991. We get global multiyear 
mean temperature T^ or T^ for each one of these satellites. 
These mean temperatures, T^ and T^ , are not adjusted for 
the coupled errors. 


In table 1, we indicate for each NOAA satellite a) the 
satellite type, b) the local time of observations, c) the time 
span of observations that is used to determine global mean 
temperature over land and ocean at two local times of the 
day, and d) the months of the overlap data that are used. In 
table 2a, we have listed for each NOAA satellite the 
unadjusted global mean Ch 2 temperatures T^ and T^, 
calculated according to equation 2 over global land. 
Similarly, in table 2b, we have listed the unadjusted global 
mean Ch 2 temperatures deduced over ocean. The differ- 
ence in the global temperature measured at two local times 
of the day, (T^-T^J, is also shown in these tables. 


Table 1. Details of satellite data used in the present analysis. 


Satellite 

Satellite 

Type 

Local Equatorial 
Crossing Times 

Time Span 
of Data 

Overlap 

Months 

NOAA 6 

Morning 

7:30 a.m. 
7:30 p.m. 

Jan. 1, 1980, to 
Dec. 31, 1981 
2 years 


NOAA 7 

Afternoon 

2:30 a.m. 
2:30 p.m. 

Jan. 1, 1982, to 
Dec. 31, 1984 
3 years 

Jan. 1982 to 
Dec. 1982 
1 2 Months 

NOAA 9 

Afternoon 

2:30 a.m. 
2:30 p.m. 

Jan. 1, 1985, to 
Dec. 31, 1986 
2 years 

Dec. 1984 and 
Jan. 1985 
2 Months 

NOAA 10 

Morning 

7:30 a.m. 
7:30 p.m. 

Jan. 1, 1987, to 
Dec. 31, 1988 
2 years 

Dec. 1986 and 
Jan. 1987 
2 Months 

NOAA 11 

Afternoon 

2:30 a.m. 
2:30 p.m. 

Jan. 1, 1989, to 
Dec. 31,1991 
3 years 

Jan. 1990 to 
Dec. 1990 
12 Months 
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Table 2a. Unadjusted average Ch 2 temperature (K) data over global land. 
(Add 200 to the temperature given in table.) 


Satellite 

No. of 
Years 

Ch 2 Average Unadjusted Temperature 
7:30 a.m. 7:30 p.m. 2:30 a.m. 2:30 p.m. 

T — T 



T 

am 

T 

PM 

T 

AM 

T 

PM 


NOAA6 

2 

51.20 

51.71 



0.51 

NOAA7 

3 



50.53 

52.04 

1.49 

NOAA9 

2 



50.41 

51.86 

1.45 

NO A A 10 

2 

51.20 

51.68 



0.48 

NOAA 11 

3 



50.54 

51.99 

1.45 


Table 2b. Unadjusted average Ch 2 temperature (K) data over global ocean. 
(Add 200 to the temperature given in table.) 


Satellite 

No. of 
Years 

Ch 2 Average Unadjusted Temperature 
7:30 a.m. 7:30 p.m. 2:30 a.m. 

2:30 p.m. 

T — T 



T 

AM 

T 

T 

* AM 

T 


NOAA 6 

2 

52.04 

52.26 



0.22 

NOAA 7 

3 



51.83 

51.86 

0.03 

NOAA 9 

2 



51.56 

51.57 

0.01 

NOAA 10 

2 

52.05 

52.32 



0.27 

NOAA 11 

3 



51.68 

51.66 

-0.02 
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From table 2b, we find the difference (T^-T^J of Ch 2 
global ocean temperature, obtained from afternoon satellites, 
has a small range, -0.02 to 0.03. Also, for the morning 
satellites, this difference has a small range, 0.22 to 0.27. 
Now, if we assume the range in (t^-T^s not due to long- 
term variations in the diurnal temperature cycle, we can 
infer noise in or T^. On this basis, we infer that such 
random noise in the 2-year or 3-year mean Ch 2 data does 
not exceed 0.03 K. Analysis of (T^-T^Jon global land, 
given in table 2a, also leads to a similar inference. 

When there is an overlap of the data from two successive 
satellites in the NOAA series, we can determine the coupled 
errors in and between those two satellites. By 
adjusting for the coupled errors in the data from 
successive satellites, we can construct a time record of the 
global a.m. temperature data from Ch 2 for 1980-1991. A 
similar global time series can be constructed for the p.m. 
data. In principle, from these time series, one could get the 
temperature change over the global ocean or global land 
separately for the a.m. and p.m. data sets. However, in order 
to enhance statistical strength in our analysis, we will 
combine the a.m. and p.m. data over the global land and the 
global ocean. 

Let us assign (t M ). to be the unadjusted average of (T AM ) i 
and , which is given by 

(t m ) =0.5x^Tam j +(t P m)J. (3) 

In this equation, the subscript i represents a number from 
1 to 5 that corresponds with one of the NOAA series of 
satellites (6, 7, 9, 10, or 11), respectively. Also, we choose 
(AT M )i to denote the difference between the average 
temperature of two overlapping successive satellites i and 
i- 1 . This can be expressed as 

|aTm j, = Ti - Ti -1 for i= 2 to5 . (4) 

With the aid of these symbolisms, we can represent the 
adjusted average temperature^), of a given satellite as 

tS j = (t m ).- X( At m ). • (5) 


In this expression, j is assigned a value from 2 to 5, which 
corresponds respectively with one of the NOAA series of 
satellites (7, 9, 10, or 1 1). Note that the uncorrected 
temperature in the series is adjusted with the sum of all of 
the preceding overlap errors. In table 3a, we are presenting 
the following global land temperature information for each 
satellite: ^(-AtJ^-AtJ,; andT*. Similar information 
for ocean is presented in table 3b. Note, from tables 3a and 
3b, the coupled error for the overlap data between NOAA 9 
and NOAA 7 is nearly the same (-0.26) over global ocean 
and land. This is because AT M between these two afternoon 
satellites is not affected by the diurnal cycle. This is not the 
case for AT m of the other three satellite overlaps. 

At this point, we may note that the MSU temperature T m of 
the global land during the 1 2-year period 1980-1991 is 
smaller than that over the ocean areas by about 0.7 K. This 
is simply due to the larger fractional area of ocean in the 
tropics as compared to that of land. 

2.4 Error analysis of MSU Ch 2 data 

In order to assess the amount of noise in the coupled errors , 
we present in figure 6a a plot of monthly mean brightness 
temperature difference vs. month, which is derived from data 
taken from NOAA 7 and NOAA 6 for 1982. The 12-month 
mean of this difference AT M for the global land is -0.015 K and 
for the global ocean is -0.304 K. The standard deviation of 
this difference over the same period is 0.084 K for global 
land and 0.034 K for global ocean. 



Figure 6a. Difference between MSU Ch 2 monthly mean 
temperature of NOAA 7 and NOAA 6 satellites, deduced 
from the overlap data for 1982. 
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Table 3a. Global land MSU Ch 2 temperature measurements and coupled errors (K). 
Included are observations before and after correction X (-AT„)for coupled errors . (Add 
200 to the mean temperature values, T„ and T* , given in table.) 




Land 


Satellite 

(tJ 

H!P!N 

X(-At.) 

(T‘) 

NOAA 6 

51.46 

N/A 

N/A 

51.46 

NOAA7 

51.28 

0.015 

0.015 

51.30 

NOAA 9 

51.14 

0.263 

0.278 

51.42 

NOAA 10 

51.44 

0.135 

0.413 

51.85 

NOAA 11 

51.26 

0.073 

0.486 

51.75 


Table 3b. Global ocean MSU Ch 2 temperature measurements and coupled errors 
(K). Included are observations before and after correction X(-AT m ) for coupled errors. 
(Add 200 to the mean temperature values, f„ and T* , given in table.) 




Ocean 


Satellite 

(t M ) 



(t:) 

NOAA 6 

52.15 

N/A 

N/A 

52.15 

NOAA 7 

51.83 

0.304 

0.304 

52.13 

NOAA 9 

51.57 



52.13 

NOAA 10 

52.18 

-0.285 

0.274 

52.45 

NOAA 11 

51.67 

0.417 

0.691 

52.36 
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A similar analysis of the coupled errors deduced from 
NOAA 1 1 and NO A A 10 overlap data is presented for 1990 
in figure 6b. The average AT M of the monthly mean bright- 
ness for 12 months is -0.073 K for global land and -0.417 K 
for global ocean; and the respective standard deviations are 
0.074 K and 0.03 K. From figures 6a and 6b, we note that 
the monthly mean brightness temperature difference for 
global land and ocean tends to have a crude mix of annual 
and semi-annual cyclical patterns. We suspect that this 
crude pattern results from changes in the diurnal cycle of 
temperature. In addition, it possibly could be due to the 
equatorial crossing of the Sun twice a year, and thereby a 
dependence on the Sun-Earth geometry during the course of 
the year. Furthermore, the pattern on the land is weighted 
more heavily by the lands in the Northern Hemisphere, 
while the pattern over the ocean is weighted by the oceans 
in the Southern Hemisphere. As a result, they tend to have 
an opposite phase. On the other hand, the 12-month mean 
value of the monthly mean brightness temperature differ- 
ence on land and ocean is strongly dependent on the 
intersatellite instrument calibration. 



Figure 6b. Same as figure 6a, except for NOAA 1 1 and 
NOAA 10 satellites, for 1990. 


The overlapping data between NOAA 9 and NOAA 7, and 
NOAA 10 and NOAA 9, covers a much shorter period (see 
table 1) and does not permit us to make annual statistical 
analyses as given above. For this reason, we will assume 
that for these two short overlap periods, the standard 
deviation, or the noise, introduced on an annual basis into 
AT M of the global land and global ocean is 0.084 K and 
0.034 K, respectively. However, since each of these short 
overlaps spans only 2 months, we will scale upwards these 
noise values by a factor of \12 / 2 . This scaling increases 


these values to 0.206 K for the global land and to 0.083 K 
for the global ocean. 

Assuming that the noise in AT M of the four overlap periods 
is random, we can calculate the effective noise in our 
satellite record for the long-term of 1980-1991. This 
effective noise in the Ch 2 temperature record calculated for 
the global ocean (ojJ is 0.063 K and for the global land is 
0.156 K. When these noise values are weighted appropri- 
ately for the land and ocean areas, we find an effective 
noise value for the globe is 0.103 K. 

As indicated earlier, the random error in the 2-year or 
3-year averages of Ch 2 temperature, deduced from each 
satellite in our analysis, is 0.03 K. In addition, based on the 
analysis of Prabhakara et al. (1995 and 1996), we assume 
the combined noise due to the hydrometeors and surface 
emissivity in these average temperatures is close to 0.05 K. 
Now, combining these errors with a ' oy c l, and yields the 
total error in the temperature derived for ocean, land, and 
the globe in the period 1980-1991. These total errors, for 
oceans o o , lands o L , and the globe ° g , can be expressed as 

O = Jo .03 2 + 0 .05 2 + 0' 2 = Vo .03 2 + 0 .05 2 + 0 .063 2 « 0 .09 K 
o \ 0 

a L = ^o”.03 2 + 0 .05 2 +C^ 2 = Vo .03 2 + 0 .05 2 + 0 .156 2 « 0 .17 K (6) 

a = ,lo .03 2 + 0 .Q5 2 + O ' 2 = Vo.03 2 + 0.05 2 + 0.103 2 « 0 .12 K . 

3 ] g 

When the satellite data is divided into a.m. and p.m. 
subsets, the corresponding errors for these subsets will be 
increased substantially. For this reason, we have not done 
the global temperature analysis for MSU Ch 2 data for a.m. 
and p.m. separately. 


2.5 Global temperature change and trend from 
MSU Ch 2: 

Now, if the initial and final temperatures of the globe are 
respectively To and Tf , then the global temperature trend 
can be deduced from the following equation: 


— r CtT „ AT ~ 

f =To + I dt = T0+ 8t , (?) 

J Hi- Al - 


where (AT/At) is the average global temperature trend over 
the period 8t. Thus, 


AT ) T f - T 0 
At j 8t 


( 8 ) 
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The adjusted global ocean temperature T „ values, shown in 
table 3b, indicates an increase of Ch 2 temperature of about 
0.21 ±0.09 K; while on the global land, this increase is 
about 0.29 + 0.17 K (table 3a). From tables 3a and 3b, we 
note that the corrected temperature over global land and 
ocean increases or decreases simultaneously during the time 
period 1980-1991. Finally, from the above analysis, we 
find that the global temperature change is 0.23 +0.12 K, 
which is based on the global mean temperature observed by 
Ch 2 (53.74 GHz) of MSU that is flown on each successive 
NOAA operational satellite from 1980 to 1991 (6, 7, 9, 10, 
and 11). 

Note the trend estimated from equation 8 is a simple 
average between the initial and final data in the time series. 
On the other hand, the trend based on a linear regression fit 
will give more importance to extreme values of the data in 
the time series. 

In figure 7, the adjusted MSU Ch 2 global temperature is 
shown as a function of time. In this figure, the solid lines 
correspond to the scale on the left ordinate and represent 
global mean temperature deduced from successive satellites 
(6, 7, 9, 10, and 1 1) covering the period 1980-1991. The 
MSU observed warming, 0.23 + 0.12 K, from 1980 to 1991 
can be readily appreciated from this figure. In figure 7, we 


have also presented the time series of the MSU Ch 2 global 
temperature anomalies, which are shown simply to illustrate 
the seasonal and interannual variations. The 12-year mean 
annual cycle is eliminated to get this anomaly series. These 
satellite-derived anomalies compare well with those 
obtained from conventional data. 

We have not considered the possible influence of the long- 
term stratospheric temperature change on that of Ch 2. The 
study of SC (1992) indicates that a decrease of about 0.3 K 
in the stratospheric temperature during the period 1979- 
1990 increases the Ch 2 temperature by only -0.02 K in that 
period. For this reason, we have neglected this small effect. 

Some NOAA satellites had drifts in their orbits, which 
change their equatorial crossing times. Intersatellite 
adjustments for coupled errors obtained from the overlap 
data of successive satellites can compensate partially for the 
errors in Ch 2 data due to such drifts. Also, the MSU 
instruments could have some degradation in their accuracy 
with time. Error due to degradation effect also can be 
remedied partially with the intersatellite adjustments for 
coupled errors. However, where there is no overlap 
between successive satellites, as is the case with the 
beginning or the ending satellite in the series, such errors 
will remain in the data. 


MSU Ch 2 Detected Global Wanning 
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Figure 7. Corrected MSU Ch 2 global temperature as a function of time. Solid lines 
correspond to the scale on the left ordinate and represent global mean temperature 
deduced from successive satellites (6, 7, 9, 10, and 11) covering the period 1980-1991. 
The MSU-observed wanning from 1980-1991 is 0.23+0.12 K (refer to text). Open 
circles correspond to the scale on the right ordinate and denote monthly temperature 
anomalies with respect to the 12-year mean annual cycle. 
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3. CONCLUSIONS 

Systematic calibration errors, present in the MSU Ch 2 data 
sets collected from sequential, polar-orbiting, Sun-synchro- 
nous NOAA operational satellites, are coupled to the 
diurnal temperature cycle over the globe. Since these 
coupled errors in MSU data differ between successive 
satellites, it is necessary to make compensatory adjustments 
to these multisatellite data sets in order to determine long- 
term global temperature change. With the aid of the 
observations during overlapping periods of successive 
satellites, we can determine such adjustments. The long- 
term time series of MSU Ch 2 global temperature adjusted 
in this fashion can yield global temperature trend. 

The error in the MSU Ch 2 global mean temperature, 
derived from each NOAA operational satellite, is -0.03 K. 
We arrive at this estimate based on the analysis of the 
difference between T m and T^ , which corresponds to 2-year 
or 3-year averages of the MSU Ch 2 temperature over the 
global ocean and land (see tables 2a and 2b). The errors in 
the intersatellite adjustments that compensate for coupled 
errors , estimated from 12 months of overlap data between 
successive satellites, are about 0.034 K for global ocean and 
0.084 K for global land areas. These errors are magnified 
when the overlap period is shorter. Such errors from short 
overlap periods, present in the Ch 2-derived temperature 


time series, effect adversely the accuracy of the estimates of 
long-term global temperature change. The combined error 
in the 1980-1991 global temperature record is -0.12 K. 

This combined error results from a) error in the Ch 2 
estimate of global mean temperature from each satellite, b) 
error introduced by hydrometeors and surface contamina- 
tion, and c) coupled errors from each overlap of successive 
satellites. 

The global temperature change is estimated from a series of 
2-year or 3-year averages of Ch 2 global temperature, which 
have been adjusted for coupled diurnal effects and 
intersatellite instrument calibration errors. We find such 
temperature change from 1980 to 1991 for the global ocean 
is 0.21 ±0.09, for the global land 0.29 + 0.17, and on the 
combined global land and ocean region 0.23+0.12. 

The advantages of the method of analysis developed here 
stem from our ability to separate the MSU data based on the 
two local equatorial crossing times of each satellite. This, 
together with the analysis of Ch 1 and Ch 2 temperature 
data (not the temperature anomalies ), allows us to under- 
stand and to determine explicitly the coupled errors on the 
global land and ocean. Also, since the estimations of 
warming, performed separately for the global ocean and 
global land, are found to be consistent with each other, 
confidence is reinforced in our method of analysis and 
its results. 
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